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ABSTRACT 



Aims. The behavior of the dark energy equation of state (EOS) is crucial in distinguishing different cosmological models. With a 
model independent approach, we constrain the possible evolution of the dark energy EOS. 

Methods. Gamma-ray bursts (GRBs) of redshifts up to z > 6 are used, in addition to type la supernovae (SNe la). We separate the 
redshifts into 4 bins and assume a constant EOS parameter for dark energy in each bin. The EOS parameters are decorrelated by 
diagonalizing the covariance matrix. And the evolution of dark energy is estimated out of the uncorrelated EOS parameters. 
Results. By including GRB luminosity data, we significantly reduce the confidence interval of the uncorrelated EOS parameter whose 
contribution mostly comes from the redshift bin of 0.5 < z < 1.8. At high redshift where we only have GRBs, the constraints on the 
dark energy EOS are still very weak. However, we can see an obvious cut at about zero in the probability plot of the EOS parameter, 
from which we can infer that the ratio of dark energy to matter most probably continues to decrease beyond redshift 1.8. We carried 
out analyses with and without including the latest BAO measurements, which themselves favor a dark energy EOS of w < - 1 . If they 
are included, the results show some evidence of an evolving dark energy EOS. If not included, however, the results are consistent with 
the cosmological constant within l<x for redshift < z < 0.5 and 2cr for 0.5 < z < 1.8. 

Key words, cosmological parameters - supernovae: general - Gamma rays: bursts 



1. Introduction 

Unexpected accelerating expansion of the universe w as first dis- 
covered by observing type l a supernovae (SNe la) (Ries s et all 
1998: iPerlmutter et al.lll999l) . This acceleration is attributed to 
dark energy, whose presence was corroborated later by other in- 
dependent sou rces including the W MAP and other o bservations 
of the CMB (ISoergel et al.ll2003l) . X-ray clusters dAllen et alj 
2002). etc. W ith more observational data available (e.g. 



Hawkins et af]|200l: lAbazaiian et al.ll200l ISpergel ei1iil [2007; 
Riess et alj 120071: IWood-Vasev et al.l l2007t iDavis et al l 2007b 
Schaeferi 120071 iPercival et al] 120071: iKomatsu et al.l [2008; 
Dunklev et ai1l2008l) . we are getting more stringent constraints 



on the nature of dark energy; nevertheless, the underlying 
physics of dark energy remains mysterious. In addition to 
the cosmological constant, many other dark energy mod- 
els have been suggested, including models of scalar fields 
(see ICopeland et al] d2006l) for a recent revie w) and modifi- 
cation of general re l ativity (see for example iDeffavetl 1 20011: 

ti ' ^ i l I In jr ^ I ■ 11 _ ^ _i I 



Binetruv et ail 120001: iMaartensI 120071: ICapozziello et all 2003 ; 
Dvali et all 120001: ICarroll et all 120041: iNoiiri & Odintsovl [2003 
2006T) . 



dard candles, have played an important role in constraining 
co smological parameters. We now have 192 samples of SN 
la dRiess et al.ll2007l: IWood-Vasev et alj|2007h iDavis et al.|[2007h 
that can be used to determine the expans ion history. And the 
proposed SNAP satellitf] (see for example Aldering et al. 2004) 
will add about 2000 samples per year. Increasing SN la sam- 
ples will provide more and more precise description of the cos- 
mic expansion. However, the redshift of the present 192 SNe 
la ranges only up to about 1.7 and the mean redshift is about 
0.5. They cannot provide any information on the cosmic ex- 
pansion beyond redshift 1.7. Here gamma-ray bursts (GRBs) 
come in and fill the void. With their higher luminosities, GRBs 
are visible across much greater distance s than supern ovae. The 
presently available 69 compiled GRBs (Schaefer 2007) extend 
the redshift to z > 6 and the mean redshift is about 2.1. After be- 
ing calibrated with luminosity relations, GRBs may be used as 
standard candles to provide information on cosmic expansion at 
high redshift and, at the same time, to tighten the cons traints on 
cosmi c e xpansion at low r edshif t. See, for example, Dai et all 
d2004l). iGhirlanda et all d2004l). iDi Girolamo etail d2005l) 



Firmani etalJ d2005l). iFriedman & Blooml d2005l). lLamb et al 
(120051). IXu et al.rd2005l). Wans & Pail 



20051) . 



Measuring the expansion history directly may be the best 
way to constrain the properties of dark energy. To measure the 
expansion history, we need standard candles at different red- 
shifts. SNe la, which are now viewed as nearly ideal stan- 



20061) . 



Su et alj d2006l) . [Schaefer "(20071), 



See http://snap.lbl.gov/ 
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IWri^htl d2007l) . and lWang etaTI ([2007) for works on GRB cos- 
mology. 

Among parameters that describe the properties of dark en- 
ergy, the equation of state (EOS) is the most important. Whether 
and how it evolves with time is crucial in distinguishing dif- 
ferent cosmological models. Due to not understanding of the 
behaviors of dark energy, simple parame tric forms such as 
w(z) = wo + w'z dC ooray & Huterer 1999 ) and vy(z) = wo + 
w a z/(l + z) d"Chevallier & Polarskil 120011: iLinder! |2003) have 



been proposed for studying the possible evolution of dark en- 
ergy. However, a simple parameterization itself greatly restricts 
the allowed wandering of w(z), a nd is equivalent t o a strong 
prior on the nature of dark energy dRiess et al.l l2007). To avoid 
any strong prior before comparing data, one can utilize an al- 
ternative approach in which uncorrected estimates are made 
of discret e w(z) of different redshifts. This approach was pro- 
posed by IHuterer & Starkmanl ((2003) and IHuterer & Cooravl 
( 200 5b and has been adopted in previou s analyses using SNe 
la dRiess et al1l2007l:ISullivan et al.ll2007al) . 

I n this work, w e apply this approach to GRB luminosity 
data dSchaeferl2007l). in addition to SN la data dRiess et alj2007t: 
IWood-Vasev et al.ll2007tlDavis et al.ll2 007). and compare our re- 
sults with those i n the previous work th at does not include GRB 
luminosity data ( Sulliv an" et alj 1200 7 a). We first briefly review 
the techniques for uncorrected estimates of dark energy evo- 
lution in section [2] The observational data and how they are in- 
cluded in the data analysis are described in section[3] We present 
our results in section|4] followed by a summary in section|5] 



2. Methodology 

Standard candles impose constraints on cosmological param- 
eters essentially through a comparison of the luminosity dis- 
tance from observation with that from theoretical models. 
Observationally, the luminosity distance is given by 



1/2 



\AttFI 



(1) 



where L and F are the luminosity of the standard candles and 
the observed flux, respectively. Theoretically, the luminosity dis- 
tance d L (z) depends on the geometry of the universe, i.e. the sign 
of Qjt, and is given by 



ddz) 



(1+z)— x 

no 



lorn ifQ * = ° 



(2) 



ll/2 



where 

E(z) = [Q,„(l + z) 3 + n x f(z) + Q t (l + z) 2 ] 1 

Q. m + Q, x + Q, k = 1 

and 



/(z) = exp 



r i + wc 

Jo 1+z 



dz 



(3) 



(4) 



Dark energy parameterization schemes enter through /(z). For 
the case where EOS is piecewise co nstant in redshift, /(z) can 
be rewritten as dSuUivan e t al. 2007a|) 



f(z„-i <z<z„) = (1 +z) 



3(l+w„) 



i=0 



3(wt— Wi+i) 



where vv,- is the EOS parameter in the / th redshift bin defined by 
an upper boundary at z,-, and the zeroth bin is defined as zp = 0. I n 
order to compare with previous analysis (Sullivan et al. 2 007al) . 
we define the first three re dshift bins to be the same as those used 
bv lSullivan etal.1 d2007ab by setting z\ = 0.2, zi = 0.5, and Z3 = 
1.8. The fourth bin is defined by z\ =7 to include GRBs. We 
carry out our analyses under two different assumptions about the 
high redshift (redshift greater than za — 7 in our case ) behavior 
of dark energy, i.e. the so-called (see lRiess et al. 1 1200% ' 'weak" 
prior, which makes no assumptions about w(z) at z > 7 and the 
"strong" prior, which assumes w(z) = — 1 at z > 7. 

In this paper we adopt x 1 statistic to estimate parameters. 
For a physical quantity £ with experimentally measured value 

standard deviation <x^, and theoretically predicted value £t(6), 
where 6 is a collection of parameters needed to calculate the the- 
oretical value, the x 2 value is given by 



cr- 



and the total y 1 is the sum of all^s, i.e. 



X 



(6) 



(7) 



The likelihood function is then proportional to exp(-;f 2 (#)/2), 
which produces the posterior probability when multiplied by the 
prior probability of 9. In the case of our analysis, the calcula- 
tion of ^- 2 s for different observational data is described in sec- 
tion|3] According to the posterior probability derived in this way, 
Markov chains are generated through the Monte-Carlo algorithm 
to study the statistical properties of the parameters. In this paper, 
we focus on the EOS parameters by marginalizing the others. 

As mentioned above, in the process of constraining cosmo- 
logical parameters, standard candles play this role by providing 
the luminosity distances at certain redshifts. However, the lumi- 
nosity distance depends on the integration of the behavior of the 
dark energy over redshift, so the estimates of the dark energy 
EOS parameters w, at high redshift depend on those at low red- 
shift. In other words, the EOS parameters vv, are correlated in the 
sense that the covariance matrix, 



C = {ww 1 ) - (w)(w T ), 



(8) 



is not diagonal. In the above equation, the w is a vector with com- 
ponents w, and the average is calculated by letting w run over the 
Markov chain. We can obtain a set of decorrelated parameters vv, 
through diagonalization of the covariance matrix by choosing an 
appropriate transformation 

w = Tw. (9) 

There can be different choic es for T. In this paper w e use the 
transformation advocated bv IHuterer & Cooravl d2005l) (see be- 
low). First we define the Fisher matrix 

F = C 1 = O t AO, (10) 
and then the transformation matrix T is given by 

T = O t AW, (11) 
except that the rows of the matrix T are normalized such that 



1. 



(12) 



The advantage of this transformation is that the weights (rows of 
T) are positive almost everywhere and localized in redshift fairly 
well, so the uncorrelated EOS param eters vv,- are easy to interpret 

' " " ~ ~^2005l) . 



intuitively dHuterer & Coora; 
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3. Observational data 

To constrain the dark energy EOS, we have made use of obser- 
vational data described below. 



3.1. Type la supernovae 

Recently compiled SN la data dRiess et all l2007t 
IWood-Vasev et all 120071; Ibavis et ail 120071) include 4 5 nearby 



Distribution of GRBs (69 GRBs in total) 



super novae dHamuv et al.l 119961 : iRiess et all 119991: JhaetalJ 



2006]), 60 ESSENCE supemovaelWood-Vasev et alJ 120071) . 
57 SN LS supernovae dAstier et al.l 120061). and 30 HST super- 
novae dRiess et al. I l200i 7 Figure Q] shows the distribution of 
these SN la samples versus redshift. The y 2 value for SNe la is 



Distribution of SNe la (1 92 SNe la in total) 
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Fig. 1. Distribution of SN la samples versus redshift 



*SN 



-z 



(Pp,i - AW) 



(13) 



where yu c , and fi p j are the observed and theoretically predicted 
distance modulus of SN la, which is defined by fi — 5 log + 25 
with the luminosity distance di in unit of megaparsec and cr lnJ is 
the intrinsic dispersion. 

3.2. Gamma-ray bursts 

Besides SNe la, GRB luminosity data is another main obser- 
vational constraint we used. As mentioned before, GRBs are 
complementar y to SNe la at high redshifts. We include GRBs 
presented by ISchaeferl d2007) (see Figure [2] for the distribution 
of these GRBs versus redshift) in our analysis by utilizing the 
five luminosity relations, i.e. the connections between measur- 
able parameters of the light curves and/or spectra and GRB lu- 
minosity: Tiag-L, V-L, E peak -L, Epe^-Ey and t rt -L 



log 
log 
log 



1 erg s 
L 

1 erg s _1 
L 

1 erg s _1 

Ey 

log^ 
1 erg 



= a\ + b\ log 
= «2 + bi log 

= A3 + £>3 log 
fl4 + Z?4 log 



Tlag(l +Z)' 



0.1 s 




v(i+ z y 




0.02 




£peak(l +Z) 


300 keV 


E ?ea kd + z)~ 



300 keV 



(14) 
(15) 
(16) 
(17) 



Maximum redshift: 6.60 
Mean redshift: 2.10 



Fig. 2. Distribution of GRB samples versus redshift 



log 



1 erg s 



a 5 + b 5 log 



trt(1 +z) 



0.1 s 



(18) 



Throughout this paper, by GRB luminosity data we refer to the 
GRBs' observational data related to such luminosity relations. It 
is worth me ntioning that the se relations may be correlated. As 
discussed in Schaefet] d2007l) . there is one significant correlation 
between the V-L and trj-L relations with the correlation coeffi- 
cient equaling 0.53. However, even for this correlation, ignoring 
it only causes a 4% unde restimate in the standard error of the av- 
erage distance modulus (Scha efer1l2007l) . so in our analysis we 
safely ignore the correlations and simply add the contributions 
from each relation (see Eq. ( l2Tb below). 

There are significant differences between SNe la and GRBs 
on the calibration. For SNe la, the calibration is done with nearby 
events and is therefore independent of cosmological parameters. 
The luminosity relations obtained in the calibration are applied 
to high-redshift events to derive the luminosity of SNe la, then 
used to constrain cosmological parameters. In this procedure, the 
calibration and the constraining of cosmological parameters are 
done separately. In contrast to SNe la, to constrain cosmolog- 
ical parameters using GRBs, we need to know the luminosity 
relations of GRBs (Eq. (1141118b ). i.e. to know the values of ct\- 
05 and b\-bs\ consequently, we need the luminosity L and the 
total collimation-corrected energy E y of GRBs, which are con- 
verted respectively from the bolometric peak flux fboio and the 
bolometric fluence Sboio of GRBs through the relations 



L 



4nd L P boio , 

Ey^oEbeam 



earn ■ 



(19) 
(20) 



The conversion depends on cosmological parameters because 
the luminosity distance di depends on cosmological models. As 
a result, the calibration and the constraining of cosmological pa- 
rameters are mixed for GRBs; i.e., we need to simultaneously fit 
calibration parameters of GRBs and cosmological parameters. 

Based on the above discussions, the x 1 value for GRBs is 
calculated by 



2 V K 



■a\ -fcilog 



TW,i(l+Zi) 



0.1 S 



■z 



lo g dnF - "2 -h log [m±f]f 
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300 keV 



with r s the comoving sound horizon at recombination and 



{i°sS-«4-Mog[%^]} 2 



{ 1o Su|f 



■a 5 -b 5 log 



trt.;(1+z;) 



0.1 s 



where L, and £ r ,■ are derived using Eq. ( fT9b and Eq. (|20l >. The 
summations run over the GRBs with the corresponding luminos- 
ity in dicator observ ed. We use the systematic errors estimated 
by (Schaefer 2007) that account for the scatter of the log-log 
plots of the luminosity versus the luminosity indicators as ov 
cr 5 in our analysis. Apparently, Xgrb * s a function of calibration 
parameters a^-a^, b\-bs and cosmological parameters that enter 
through the luminosity distance d^. 

3.3. Other data 

In addition to SNe la and GRBs, we have also used the con- 
straints below follow ing previous analyses dRiess et al. I l2007t 
ISullivan et al.1l2007ah 

- Constraints on dimensionless mass densities: The SDSS 
large-scale structure measurements give the constraint 
on loc al mass density in terms of £l„,h = 0.213 ± 
0.023 dTegmarketal .1120041) . The WMAP three-year data 
combined with the HST key project c onstraint on th e Hubble 
constant gives Q. k = -0.014 ± 0.017 dSpergel et al.ll2007l) . 

- The SDSS luminous red galaxy, baryon acoustic oscilla- 
tion (BAO) distance parameter to Zbao = 0.35: A = 

W- PfeBAo)^^] 173 , w here r(z) = d L (z) /(l + z).A = 

0.469 (0%) 35 ± 0.017 from lEisenstein et al.1 d2005l) and 
the th ree-year WMAP results give n = 0.95 dSpergel et al.1 
120071) . 

- The distance to last scattering, z—1089: If nonzero cosmic 
curvature is a llowed as we do in o ur analysis, the three-year 
WMAP data dSpergel et al.ll2007l) gives the shift parameter 

Rcmk = ^Pmzcmb) = 1.71 ± 0.03 dWang & Mukheried 
120071) . 

- The distance ratio between zbao = 0.35 and zcmb = 1089: 



[r 2 (ZBA0)«^] 



1/3 



-Ko.35 = : : ■ (22) 

Kzcmb) 

The SDSS BAO analysis dEisenstein et al.l 120051) gives 
R . 3 5 = 0.0979 ± 0.0036. 

The corresponding x 2 f° r these constraints are directly calcu- 
lated using Eq. (0. 

We have also studied the dark energy EOS evolution with 
the above BAO con straints replaced by the latest BAO measure- 
m ents presented in iPercival et al.l d2007l) . for which the x 1 value 
is dPercival et al.ll2007l) 



X\,KQ - ^BAO^BAO^BAO 

where 



D v (0.35) 



0.1094 J 



(23) 



(24) 



1 1 35059 -24031 
c bao - I _24031 108300 /' 



(25) 



Thi s constraint itself favors a dark energy EOS of w < 
-1 dPercival et al.ll2007l) . 



(21) 4. Results 



Figures [3] and [4] show our results for the weak prior and 
strong prior respectively. For these two figures, we have in- 
cl uded subsets o f data f rom section [331 same as that are used 
in ISullivan et al.l (l2007al) besides SNe la. For the results pre- 
sented in Figure [5] t he BAO constraints are updated with the lat- 
est measurements dPercival et alj |2007). see Eq. d23l , Eq. d24l . 
and Eq. (1251 . A comparison between Figures[3]and[4]shows that 
the results are insensitive to the priors, i.e. insensitive to whether 
w(z > 7) = -1 is assumed or not for dark energy. 

Since Figures [3] and [4] only differ from results derived 
by ISullivan et ail l2007al) in that we include GRB lumi- 
n osity data, compariso ns of Figures [3] and H] with figures 
in lSullivan et alj (12007a) demonstrate the improvement made by 
including GRBs. We find that there is little improvement in w\ 
and W2- This is because at low redshift, where we have both SNe 
la and GRBs, there are fewer GRBs than that of SNe la (see 
Table [T] in the first two bins the number of GRBs is negligible 
compared with that of SNe la); at the same time, the contribu- 
tions to w\ and W2 from high redshift, where we have a consid- 
erable number of GRB samples (see Table [TJ, are too small (see 
the weight histograms in Figures[3]and|4|i to improve constraints 
on w\ and W2 significantly. The most significant improvement 



bin 


1 


2 


3 


4 


redshift range 


0-0.2 


0.2-0.5 


0.5-1.8 


1.8-7 


number of SNe la 


47 


59 


86 





number of GRBs 


1 


3 


32 


33 


total number 


48 


62 


118 


33 



Table 1. Number of SNe la and GRBs that fall into the four bins 



lies in W3, whose contribution mostly comes from the third bin, 
where we have several GRBs (see Table [1]). The lcr confidence 
interval of W 3 with GRBs included is less than one third of that 
presented in Sulliv an et al.l (l2007al) without including GRB lu- 
minosity data. 

For Figures [3] and |U w\ and wi are consistent with the cos- 
mological constant within lcr, and W3 consistent within 2cr. 
While in Figure [5] for which the latest BAO measurements are 
used instead, the cosmological constant lies outside of the lcr 
confidence intervals of w\ and W2, and outside the 2cr confidence 
interval of W3, though still inside the 2cr confidence intervals of 
wi and W2- These results show some evidence of an evolving 
dark energy EOS. This is not surprising provided that the lat- 
est BAO measurements thems elves favor a dark energy EOS of 
w < -1 dPercival et al.l 120071) . The BAO distance information 
lies in the second redshift bin, so including it leads to a smaller 
h>2. And main data we used depends on the integration of the 
dark energy evolution, thus the decrease in W2 causes increases 
in vvi and W3 . 

The constraints on w\ are very weak. The uncertainty is so 
great that we plot its probability separately. This is due to three 
reasons. First, there are not enough samples of standard candles 
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Fig. 3. Estimates of the uncorrelated dark energy EOS parame- Fig. 4. Estimates of the uncorrelated dark energy EOS parame- 
ters using the weak prior. In turn are the plots of vv, (z = 1-3) ters using the strong prior. Same as Figure [3] except using the 
versus redshift, window functions of vv, (z = 1-4) with respect strong prior, 
to the 4 bins, probability distribution of vv,- (z = 1-3), W4, and 
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Fig. 5. Estimates of the uncorrelated dark energy EOS pa- 
rameters using the strong prior. Same as Figure |4] ex- 
cept BA O constraints wer e updated with the latest measure- 
ments dPercival et alj|2007l) . 



in the fourth bin, all of which are GRBs. From Table Q] it can 
be seen that the number ratio of third bin to the fourth bin is 
about 4. Second, as mentioned earlier, the estimate of the be- 
havior of dark energy at high redshift depends on its behavior at 
low redshift; consequently, the uncertainty of EOS parameters 
at low redshift will be reflected on EOS parameters at high red- 
shift. Therefore we get increasing errors as the redshift increases. 
Thirdly, the density ratio of dark energy to matter is given by (as- 
suming a constant EOS parameter for dark energy) 

Px_ _P*o(l+z) 3(1+w * ) 

pm 



3(1 +zf 



(26) 



Pm0(l +Z) 3 

For negative w x , the ratio decreases as z increases. For example, 
when w x = - 1 , the ratio is about 1 /9 at z — 2 . At higher redshift, 
matter dominates over dark energy, then dark energy becomes 
less important in determining the cosmic expansion. Thus the 
constraints imposed on the behavior of dark energy by the ex- 
pansion history become weak compared with that at low redshift 
where dark energy is important. Despite the large uncertainty in 
H>4, there is indeed some restriction imposed by GRBs. From the 
probability plots of w\ in Figures [3] S and[5] it can be seen that 
there is obviously a cut at about zero. In other words, it is most 
probable that the ratio in Eq. ( |26| i continues to decrease at a red- 
shift beyond 1 .8. The probability cut at the left of -200 is due to 
the precision of the computer and can be viewed as the negative 
infinity. To get substantial constraints on the dark energy EOS 
beyond 1.8, we need more GRB samples. 

To see the overall improvement made by including GRB lu- 
minosity d ata, we calculate the fi gure of merit (FOM), which is 
defined by dSullivan et al.ll2007allbl) 



FOM = 



nl/2 



(27) 



For the the results presented in Figure|5] FOM = 9.6. And if the 
GRB luminosity data are excluded, FOM = 8.8. 

5. Summary 

We used a model-independent approach to constrain the evo- 
lution of dark energy. First, we separated the redshifts into 4 
bins and assumed a constant EOS parameter for dark energy 
in each bin, then estimated the uncorrelated EOS parameters. 
We mainly used the SNe la and GRBs in our analysis. Other 
constraints from SDSS, 2dFGRS, HST, and WMAP are also in- 
cluded. Compared with the results obtained without including 
GRB luminosity data, the confidence interval of the third uncor- 
related EOS parameter, whose contribution mostly comes from 
the third bin, is reduced significantly. Even though constraints 
at high redshift where we have only GRBs are very weak, from 
the obvious probability cut of the EOS parameter at about zero, 
we can infer that it is most probable that the ratio of dark energy 
to matter continues to decrease beyond redshift 1.8. To get sub- 
stantial constraints at redshifts beyond SNe la more GRBs are 
needed. 

If the latest BAO measurements, which themselves favor a 
dark energy EOS of w < -1, are included, the results show 
some evidence for an evolving dark energy EOS. Otherwise, the 
results are consistent with the cosmological constant. 
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